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Summary  of  Work 

The  work  accomplished  to  date  can  be  divided  into  two  categories: 

(1)  component  models  and  (2)  model  integration.  The  latter  is 
accomplished  with  a general  purpose  process  simulator  program  for 
water  reuse  processes. 

Component  Models 

The  component  modules  include  the  following:  (1)  ultrafiltration, 

(2)  reverse  osmosis,  (3)  ozonation,  and  (4)  hypochlorination.  The 
simulator  structure  enables  a model  for  the  equalization  tank  to 
be  created  using  the  general  features  of  the  simulator. 

Ultrafiltration 

The  present  model  for  the  tubular  ultrafiltration  module  is  based 
on  the  film  theory  for  mass  transfer  through  the  boundary  layer 
at  the  inner  wall  of  the  UF  tube.  Permeability  coefficients  are 
included  for  water,  dissolved  solids,  and  TOC.  For  a given  module 
configuration,  inlet  flow,  inlet  suspended  solids  concentration, 
inlet  dissolved  solids  concentration,  inlet  TOC  concentration,  and 
pressure  drop  across  the  module,  the  model  predicts  the  permeate 
flow,  permeate  dissolved  solids  concentration,  and  permeate  TOC 
concentration.  The  model  assumes  that  no  suspended  solids  pass 
through  the  membrane.  This  model  is  described  in  detail  in  the 
interim  report  "A  Mathematical  Model  of  a Tubular  Ultrafiltration 
Unit  for  Water  Re-Use  Systems". 

Reverse  Osmosis 

The  model  for  the  hollow  fiber  reverse  osmosis  module  is  also  based 
on  the  film  theory  for  mass  transfer  through  a boundary  layer  formed 
at  the  outer  surface  of  the  hollow  fibers.  For  a given  module 
configuration,  inlet  flow  rate,  inlet  dissolved  solids  concentration, 
inlet  TOC  concentration,  and  pressure  drop  across  the  module,  the 
permeate  flow,  permeate  dissolved  solids  concentration,  and  permeate 
TOC  concentration  are  predicted.  Analogous  values  for  the  concen- 
trate stream  are  calculated  by  material  balance.  This  model  is 
described  in  detail  in  the  interim  report  "A  Mathematical  Model  of 
a Hollow  Fiber  Reverse  Osmosis  Unit  for  Water  Re-Use  Systems". 

Ozonation 

The  model  for  ozonation  in  the  presence  of  UV  light  is  based  on 
basic  kinetic  mechanisms  for  the  reaction  and  decomposition  of  ozone. 
The  ozone  enters  as  a gas,  is  dissolved  Into  the  aqueous  phase  and 
then  either  reacts  with  the  organic  material  in  the  water  or  decom- 
poses. For  a given  contractor  configuration  and  inlet  water  con- 
ditions, the  ' del  predicts  the  effluent  water  TOC  concentration. 

The  model  is  aescribed  in  detail  in  the  interim  report  "A  Mathema- 
tical Model  of  an  Ozonation  Contacting  Unit  for  Water  Re-Use  Systems" 


Hypochlorination 


The  model  for  the  hypochlorination  module  is  based  on  the  dis- 
sociation mechanism  for  hypochlorite.  The  model  predicts  the 
free  and  available  chlorine  concentration  from  inlet  stream  con- 
ditions. The  model  is  described  in  detail  in  the  interim  report 
"A  Mathematical  Model  of  a Hypochlorination  Unit  for  a Water  Re- 
Use  System". 


Model  Integration 


The  component  modules  were  incorporated  into  the  process  simulator 
to  facilitate  the  simulation  of  a variety  of  configurations  for 
water  reuse  processes. 


amic  Process  Simulator 


In  addition  to  the  modules  specific  to  water  re-use  processes, 
the  simulator  includes  modules  for  general  process  elments  such 
as  mixed  tanks,  pumps,  stream  splitters,  etc.  The  simulator 
permits  the  configuration  of  the  water  reuse  process  to  be 
specified  via  data  cards,  along  with  specifications  for  the  input 
streams.  The  simulator  calculates  the  flow  and  concentration 
of  each  stream  in  the  process  as  a function  of  time.  The  results 
may  be  presented  in  tabular  form  or  may  be  plotted.  To  provide 
some  assurance  that  the  models  and  process  configuration  is  con- 
sistent, total  and  component  material  balances  are  computed  for  the 
run.  The  dynamic  process  simulator  is  described  in  detail  in  the 
interim  report  "A  Dynamic  Model  for  Application  in  Water  Re-Use 
Systems". 


Integrated  Model  Example 


Using  the  dynamic  process  simulator,  a simulation  run  was  made 
using  the  configuration  proposed  for  the  MUST  WPE.  The  run  demon- 
strated that  the  component  models  behaved  as  expected  in  an  inte- 
grated configuration,  and  that  the  dynamic  process  simulator  could 
successfully  simulate  a realistic  water  re-use  process.  Simulation 
of  a 20-hour  cycle  required  24  minutes  of  computer  time  on  an  IBM 
370/3033. 
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2.0  Publications  and  Presentations 


The  following  publications  and  presentations  have  resulted  from 

this  work: 

1.  Starks,  D.  M. , and  C.  L.  Smith,  "A  Dynamic  Model  of  a Water 
Treatment  Unit",  Proceedings  of  the  1977  IEEE  Conference  on 
Decision  and  Control,  New  Orleans,  December  7-9,  1977. 

2.  Starks,  D.M. , and  C.  L.  Smith,  "A  Mathematical  Model  of  a 
Reverse  Osmosis  System  for  Water  Purification  Use",  Proceedings 
of  the  Ninth  Annual  Pittsburgh  Conference  on  Modeling  and 
Simulation,  Pittsburgh,  April  27-28,  1978. 

3.  Starks,  D.M. , and  C.  L.  Smith,  "A  Mathematical  Model  of  an 
Ultrafiltration  System  for  Water  Purification  Use",  Proceedings 
of  the  1978  Summer  Computer  Simulation  Conference,  Newport 
Beach,  July  24-26,  1978. 

Copies  of  all  of  these  papers  are  included  in  the  appendix. 
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Abstract 

For  support  of  field  medical  units,  the  U.S. 
ray  Is  developing  a self-contained  waste  water 
treatment  system  to  produce  potable  water  for  use 
within  the  unit.  To  aid  In  developing  the  control 
strategy  and  the  fault  detection/fault  Isolation 
logic,  a dynamic  model  of  the  system  Is  being  de- 
veloped. 

The  first  phase  of  the  effort  Is  to  develop 
models  of  the  system's  conponents,  such  a ultrafil- 
tration, reverse  osmosis,  ozonation,  etc.  These 
models  are  then  combined  to  obtain  a model  for  the 
complete  system.  To  the  greatest  extent  possible, 
fundamental  equations  form  the  basis  of  the  models. 

In  determining  coefficients  for  the  models,  the 
the  data  was  obtained  from  the  dxperlments  that 
formed  the  basis  of  the  design  of  a pilot  plant  ver- 
sion of  the  unit.  To  refine  the  model,  this  data 
Is  chen  supplemented  by  operational  data  from  the 
pilot  plant. 


1.  Introduction 

The  U.S.  Army  has  a requirement  to  provide  a 
mission-oriented  medical  treatment  system  which  Is 
designed  and  equipped  to  facilitate  rapid  estab- 
llahmenc  and  disestablishment.  The  flexibility 
permits  ltaedlate  response  for  a medical  support 
unit  to  any  tactical,  environmental  or  geographical 
change.  This  system  will  provide  a contamination- 
free  and  controlled  environment  in  which  medical 
surgical  and  ancillary  procedures  and  other  sup- 
ptrging  functions  can  be  performed.  The  moHlle 
medical  treatment  system  Is  termed  the  MUST  Medical 
Complex,  Associated  with  the  MUST  Medical  Complex 
is  a Water  and  Waste  Management  Subsystem  (WWMS). 
This  subsystem  is  required  to  treat  and  dispose  of, 
without  degradation  of  the  environment  or  danger  to 
personal  health,  all  toxic  and  contaminated  waste 
materials  generated  within  the  functional  areas  of 
the  Medical  Complex.  In  addition  to  the  waste 
treatment  and  disposal,  the  WPE  within  the  WWMS 
must  be  capable  of  producing  potable  water  from  a 
fresh  or  brackish  water  source  and  nnncomsumptlve 
reuse  water  from  the  MUST  Medical  Complex  waste 
water  effluent. 

The  objective  of  the  present  program  la  to: 

1.  Develop*  an  Integrated  dynamic  model  de- 
scribing the  operational  characteristics 
of  the  mater  prdceaslng  element.  Emphasis 
will  be  placed  on  the  reuse  mode  of  opera- 
tion utilizing  the  MUST  hospital  coaposlte 
waste  or  the  x-ray,  laboratory,  and  OR 


composite  waste.  The  methodology  to  apply 
the  model  to  other  configurations  and 
other  wastes  will  be  developed  and  doc- 
umented. 

2.  Using  the  dynamic  model  as  the  basis,  de- 
velop a control/monitoring  system  for  the 
operation  of  the  WPE. 

3.  Using  the  dynamic  model  as  the  basis,  de- 
velop a fault  detection/fault  Isolation 
package  for  the  WPE. 

This  presentation  will  only  consider  the  first  of 
these  objectives. 


2.  System  Description 

When  in  the  re-use  mode  of  operation,  the  WPE 
utilizes  the  following  components  in  series: 

1.  Equllzation  (mixing)  tank 

2.  Ultrafiltration 

3.  Microfiltration 

4.  Reverse  osmosis 

5.  Ozonation  (accompanied  by  exposure  to 
ultraviolet  light) 

6.  Clorinatlon 

Dynamic  models  are  required  for  each  of  these  com- 
ponents. 

The  sources  of  the  recycle  water  include 
showers,  operating  room,  kitchen,  laboratory,  and 
x-ray.  Due  to  the  large  number  of  sources  of  con- 
taminants, It  Is  inpractlcal  to  Identify  every 
molecular  compound.  In  such  cases,  the  usual  prac- 
tice la  to  Introduce  pseudo-counpounds,  which  will 
consist  of  several  molecular  compounds  of  similar 
characteristics.  Currently,  only  two  pseudo-com- 
pounds see  being  used,  with  all  contaminants 
being  treated  as  either  dissolved  or  suspended 
solids. 

In  developing  the  design  for  the  WPE,  the 
first  step  was  to  develop  and  operate  a 1/4-scale 
breadboard  system.  The  next  step  was  to  construct 
a full  scale  pilot  plant.  Data  from  both  serve  as 
the  experimental  basis  for  the  dynamic  rodel. 


3.  Model  Formulation 

The  dynamic  models  for  each  component  of  the 
WPE  are  based  on  fundamental  concepta,  such  aa 
heat  balancea,  material  balancea,  meat  transfer 
equatlone,  etc.  For  exapple,  the  model  for  the 
ultrafiltration  unit  involves  the  following  equa- 
tions: 
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4.  Model  Calibration 


In  order  for  the  model's  behavior  to  match  that 
of  the  physical  system,  the  unknown  parameters  In 
the  various  relationships  are  adjusted  until  the 
model  results  agree  with  the  experimental  results 
to  within  an  acceptable  tolerance.  In  the  ultra- 
flltratlon  model,  five  parameters  must  be  deter- 
mined in  this  manner.  In  this  case,  experimental 
data  Is  currently  available  only  for  steady  opera- 
tion, although  future  plans  call  for  data  to  be 
collected  during  periods  of  non-steady  state  op- 
eration. 

After  the  parameters  have  been  fit  to  the  ex- 
perimental data,  the  model  behavior  must  be  ob- 
served under  several  different  conditions  to 
verify  that  Its  behavior  appears  to  be  "reasonable". 
Futher  tests  of  this  type  are  necessary  after  the 
component  models  have  been  Integrated  Into  a model 
of  the  complete  system. 

The  presentation  will  describe  the  experiences 
of  the  authors  In  developing  component  models  for 
the  WPE,  wiht  emphasis  being  on  the  process  of  es- 
tablishing a valid  model  from  the  available  experi- 
mental data. 


the  boundary  layer  at  the  tube  wall. 

Models  of  other  elements  of  the  WPE  Involve  rela- 
tionships such  as  heat  transfer,  reaction  kinetics, 
etc. 

Although  the  basic  theoretical  equations  are 
relatively-  easy  to  develop,  numerical  values  for 
the  various  coefficients  are  more  difficult  to  ob- 
tain. Empirical  correlations  are  also  often  nec- 
essary. In  the  ultrafiltration  model,  the  follow- 
ing are  used: 

1.  The  mass  transfer  coefficient  is  related 
to  the  diffuslvity. 

2.  The  osmotic  pressure  Is  related  to  the 
concentration  of  contaminant. 

Furthermore,  some  physical  properties  (such  as  the 
diffuslvity)  may  be  functions  of  operating  condi- 
tions (such  as  contaminant  concentrations). 
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A MATHEMATICAL  MODEL  OF  A REVERSE  OSMOSIS 
SYSTEM  FOR  WATER  PURIFICATION  USE 


David  M.  Starks  and  Cecil  L.  Smith 
Department  of  Chemical  Engineering 
Louisiana  State  niversity 
Baton  Rouge,  Louisiana  70803 


ABSTRACT 

This  paper  describes  the  development  of  a mathematical  model  depicting  the  reverse  osrosis  ur. 
employed  in  the  water  processing  element  under  development  by  the  U.S.  Army  Medical  Research 
and  Development  Command  for  support  of'field  medical  units.  Methodology  for  solving  the  mode 
equations  and  fitting  the  model  parameters  to  experimental  operating  data  has  been  developed 
and  programs  written  for  a digital  computer.  This  article  presents  in  detail  the  model  eoua- 
tions,  the  parameter  fitting  procedures,  and  some  typical  results. 


DEVELOPMENT  OF  BASIC  EQUATIONS  AND  METHODS  USED 


The  model  equations  developed  characterize  the  operation  of  the  duPont  B-10  Reverse  Osmosis 
Separator  as  illustrated  in  Figure  1.  Contaminated  water  enters  via  the  high  pressure  feed 
tube.  Flowing  in  the  radial  direction,  contaminated  water  becomes  progressively  more  concen- 
trated as  water  preferentially  permeates  through  the  membrane  wall  into  the  inner  bore  of  tht 
fibers.  Concentrate  is  collected  at  the  flow  screen  and  moves  axially  towards  the  ccncentrai 
outlet.  Purified  water  collected  inside  the  fibers  exits  the  module  as  permeate. 


Figure  1.  Reverse  Osmosis  Separator 


FORMULATION  OF  EQUATIONS 

...To  develop  the  area  available  for  the  flux  of  water  and  contaminants  into  the  fiber  bore,  le 
Nf  be  the  number  of  perfectly  straight,  evenly  distributed  fibers  parallel  to  the  axial  dire 

tion  per  unit  of  cross  sectional  area.  Collectively,  these  idealized  fibers  would  have  the 
same  fluxes  of  water  and  contaminants  as  the  actual  fiber  configuration. 
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H>r  the  oil terential  element  in  figure  2.  the  area  available  for  flux  becomes: 

(Nf)(nDFL)[7i(r  + A r)2  „ *(r)2) 
where:  Df  is  the  outside  diameter  of  the  fiber  (m) 


L is  the  length  of  the  fibers  (m) 
r is  the  radial  coordinate  (m) 

A steady  state  material  balance  on  the  dissolved  solids,  component  ‘A1,  for  the  slice  of 
module  between  r and  r + Ar  in  Figure  2 yields: 


Input  rate  of  A = Output  rate  of  A 


(27irL)V(r)CA(r)  = 2m(r  + Ar)LV(r  + Ar)CA(r  + Ar)  + Jft( r ) ( Nf nDj-L ) [u ( r + Ar)2  -nr2 


where:  V(rl  = radial  velocity  of  fluid, (m/hr) 
C«(r)  = concentration  of  A (gm/m^) 


0A(r)  = flux  of  A into  inner  channel  of  hollow  fibers  (gm/hr-m; 


A similar  steady  state  material  balance  on  the  TOC,  component  'C',  for  the  differential 
element  in  Figure  2yields: 


(2arL)V(r)Cc(r)  = 2*(r  + Ar)LV(r  + Ar)Cc(r  + Ar)  + Jj-(r)  (NpirDpL)  [tt(  r + Ar)2-nr2 
where:  Cj.(r)  = Total  Organic  Carbon  (TOC)  concentration  (gm/m^) 

0c(r)  = °f  C into  inner  channel  of  hollow  fibers  (gm/hr-m2) 

A steady  state  total  material  balance  yields: 


(2itrL)V(r)p(r)  = [2n(r  + Ar)LjV(r  + Ar)p(r  + Ar)  + [ JA( r ) + Jg(r) 


+ Jc(r)](NfmDpL)[Ti(r  + Ar)2  - ur2] 


where:  p(r)  = total  mass  density  (gm/ni  ) 


JB(r)  = flux  of  water,  component  B,  into  the  fiber  bore  (gm/m  -hr) 


In  addition,  a steady  state  material  balance  on  component  ' B ' can  be  written  to  obtain  an 
equation  analogous  to  equation  (1).  However,  this  equation  is  not  necessary  as  it  become 
dependent  equation  when  equations  (1),  (2),  and  (3)  are  used  in  the  model. 


V(r+Ar) ,CA(r+Ar) ,Cp(r+Ar) 


Cc(r+Ar),p(r+Ar) 


V(r), CA(r),CB(r),p(r} 


tv-'-.! 


\ \ 


Figure  2.  Radial  Differential  Element 


ANALYTICAL  EQUATIONS. 

Equations  (1),  (2)  and  (3)  are  readily  converted  to  differential  equations  by  taking  the 
limit  as  Ar  approaches  zero.  The  results  are: 


i [rV(r)CA(r')J  = - NfliDprJA(r) 

■jj~  lrV(r)Cc(r)]  = - NfnDFrJc(r) 

[rV(r)p(r)]  = - Nf nDFrlOA(r)  + Jg(r)  + Oc(r))  { 

The  boundary  conditions  are  as  follows: 

VOV  = ^/(^RiD 
W = CAo 
CC(IV  “ CCo 

p(R^)  = calculated  from  CA(R^)  and  CC(R.) 

3 

where:  F = feed  rateto  module,  m /hr 

0 3 

C.  = inlet  concentration  of  A,  gm/m 
Ao  3 

r = inlet  TOC  concentration,  gm/m 
Co 

FLUXES  AT  THE  BOUNDARY  LAYER 

To  obtain  the  fluxps  JA,  Jg,  and  J^.,  a stagnant  boundary  layer  is  imagined  to  exist  near 

outer  surface  of  the  hollow  fiber  as  shown  in  Figure  3.  The  following  fluxes  into  and  ou 
the  boundary  layer  are  defined:’ 


■ B<CA2  - CA3> 

(7) 

JDA  = kA(CA2  •'  CA15 

(10) 

= y(AP  - An) 

(8) 

JDB  = kB^CBl  ‘ CB2^ 

(11) 

* C(CC2  - CC3) 

(9) 

JDC  = kC*CC2  ' CCp 

(12) 

where:  0pA  = diffusional  flux  of  contaminant  out  of  the  boundary  layer,  gm/m-hr 
J„D  = diffusional  flux  of  water  into  the  boundary  layer,  gm/m2-hr 

Uu  « 

Jnr  = diffusional  flux  of  TOC  out  of  the  boundary  layer,  gm/m-hr 

. 3 

cAr  = concentration  of  A in  bulk  fluid,  gm/m 

CA2  = concentration  of  A in  boundary  layer,  gm/m2 

CA 3 c concentration  of  A in  permeate,  gm/m2 

3 

C„,  * concentration  of  B in  bulk  fluid,  gm/m 

Dl  « 3 

C„„  = concentration  of  B in  boundary  layer,  gm/m 
oi  3 

C„,  = concentration  of  B in  permeate,  gm/m 

Cr , * concentration  of  C in  bulk  fluid,  gm/m 
u 3 »• 

C_0  = concentration  of  C in  boundary  layer,  gm/m 

Lc  3 '■ 

= concentration  of  C in  permeate,  gm/m 
to  2 

Y = permeability  coefficient  for  pure  water,  gm/m  -hr-atm 

B * permeability  coefficient  for  dissolved  solids,  m/hr 

C * permeability  coefficient  for  TOC,  m/hr 

AP  = pressure  drop  across  membrane,  atm 

Am  = osmotic  pressure  difference  across  membrane  » - iij,  atm 

Tig  = osmotic  pressure  at  boundary  layer,  atm 
ti3  « osmotic  pressure  of  permeate,  atm 


kA,kr  • mass  transfer  coefficients  for  solids  and  TOC,  respectively,  m/hr 
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CA3’CB3,CC3 


Figure  3.  Fluxes  into  and  out  of  boundary  layer 


The  osmotic  pressure  is  related  to  the  concentration  of  contaminant,  CE,  and  the  absolute 
temperature,  T,  by  the  empirical  relationship 

* = f(CE,T)  = aTCE(l  + 0C£)2  (13) 

where:  a = coefficient,  atm-m3/gm-°i< 

. 3 

B = coefficient,  m /gm 
T = absolute  temperature,  degrees  Kelvin 

3 

CE  = effective  contaminant  concentration  = + •tEcc»  gm/«n 

kE  = coefficient  expressing  the  concentration  of  A equivalent  to  a unit  concentrati 
L of  C for  purposes  of  osmotic  pressure,  dimensionless 

To  obtain  A*  * n2  - r3,  it  is  necessary  to  evaluate  h2  using  and  CE2;  using  C A3  and 
CC3' 

The  flux  Jp  due  to  the  bulk  flow  toward  the  boundary  layer  can  be  obtained  by  writing  a to 
material  balance  around  the  boundary  layer: 

JF  = JA  + JB  + JC  + JDA  ' JDB  + JDC  O4) 

In  addition,  compone  it  balances  for  component  A and  for  component  C can  be  written: 

JFXA1  " °A  + JDA  (15)  = Jr  + Jnr  («) 


where  Xftl,  XC1  = mass  fractions  of  dissolved  solids  and  TOC  in  bulk  stream,  respectively. 

A component  balance  can  also  be  written  for  component  B,  but  it  will  not  be  an  independent 
■equation. 

In  addition,  the  mass  fractions  X^3  and  X^3  are  related  to  Jg,  and  Jj.  as  follows: 


A3  JA+JB  + JC 


C3  JA  + °B  + °C 


Correlations  for  the  mass  transfer  coefficients,  k.  and  k.,  have  been  developed  analogous 
to  those  available  for  flow  across  heat  exchanger  *tube  bundles.  The  final  result  is: 
e,  D V 


k - C,  (N-  ) \ 


(19)  fin.  - 


nr*  W-*  UTF  ¥ * IT".1*  v ' 
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coefficient  whicn  is  a luninon  ui 
various  physical  constants 
exponent  of  the  Reynolds  number 


where:  Nge  = Reynolds  number 

0 = outer  diameter  of  fiber 

o 

V = radial  velocity 

v = kinematic  viscosity 

k = mass  transfer  coefficient 

The  same  values  for  the  exponent  ej  will  be  used  for  both  the  solids  and  TOC  but  different 
values  will  be  used  for  C^,  owing  to  different  physical  properties  of  the  two  components. 

APPROXIMATIONS  IN  THE  FLUX  EQUATIONS 


& 

§ 

p. 

S 

* 


c< 

V, 

£ 


JL 


In  the  reverse  osmosis  unit,  the  concentrations  of  dissolved  solids  and  TOC  will  be  relativ 
small.  Furthermore,  the  RO  membrane  will  tend  to  reject  these  components.  For  these  reasc 
the  following  assumptions' will  be  made: 

1.  The  membrane  permeability  fluxes  and  are  small  compared  to  Jg. 

2.  The  diffusional  fluxes  JDA>  J^g,  and  JD(-  are  small  compared  to  0g. 

3.  The  density  is  that  of  water  namely  p. 

Table  1‘  gives  the  flux  equations  with  these  approximations  included. 

Table  I.  Simplified  Equations  for  Fluxes  at  the  Boundary  Layer 

Osmotic  Pressure  Relationships 


JA  = BlCA2  ' CA3) 

Ait  = itg  * 'n3 

Og  = y(AP  - An) 

*2  = 

JC  = C(CC2  ” CC3) 

*3  = “TCE3(1 

JPA  ° kA*CA2  " CAp 

CE2  = CA2  + kl 

JDC  = kC(CC2  ' CC1^ 

CE3  " CA3  + kl 

:12)22 

cE3)2 


Material  Balance  at  Boundary  Layer 


E C2 
ECC3 


Mass  Fraction/Density/Concentration  Relationshi 


B 

JA  + JDA 


Vai 

JFXC1  " JC  + JDC 


A3 


X = — 
*C3  JD 


B2  ‘ 

1 ' XA2  ‘ XC2 

XB3  = 

1 -XA3 

A2 

PXA2 

CA3  = 

PXA3 

B2  = 

PXB2 

CB3  “ 

PXB3 

C2  * 

PXC2 

CC3  “ 

PXC3 

- X 


C3 


SOLUTION  METHOD 

The  solution  method  basically  consists  of  numerically  integrating  the  differential  equati 
resulting  from  the  total  and  component  material  balances.  However,  this  is  somewhat  comp 
cated  by  the  fact  that  the  fluxes  must  be  determined  at  each  increment  by  solving  the  set 
nonlinear  equations  summarized  in  Table  1. 

PARAMETER  ESTIMATION  CONSIDERATIONS 

Values  of  several  key  variables,  such  as  mass  transfer  coefficients  and  osmotic  pressure 
relation  coefficients,  must  be  determined  from  experimental  data.  A Pattern  Search  strat 
was  used  to  find  a set  of  values  for  these  parameters  which  minimized  a cost  function  exp 
sive  of  the  difference  between  model  calculations  and  experimental  results. 

Table  II  lists  the  parameters  whose  values  have  been  obtained  by  the  parameter  cstimatic 
procedure.  The  parameters  Nf,  y,  B.  and  C may  only  be  obtained  from  the  operation  of  the 
module  itself. 


Of  the  remaining  parameters,  three  of  these  (CJA, 


and 


ej)  are  required  for  the  mass 


transfer  relationships  and  three  (a,  6,  and  kr)  are  required  for  the  osmotic  pressure  rel 
ship.  It  should  be  possible  to  develop  experimental  procedures  to  evaluate  these  indeper 
of  the  RO  module  test,  but  in  the  absence  of  these  independent  experiments,  these  six  P*> 
meters  are  also  evaluated  from  the  RO  test  data. 
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Tabic  11  Model  Parameters 


WU>i>  J\  wV.V.V  Vc  VY  VI 33 S21 


Number  of  ideal  fibers  per  unit  cross-sectional  area  of  the  fiber  bundle 

Permeability  coefficient  for  pure  water 
Permeability  coefficient  for  dissolved  solids 
Permeability  coefficient  for  TOC 
Exponent  of  Reynolds  number  in  transfer  correlation 

Coefficient  in  osmotic  pressure  correlation 
Coefficient  in  osmotic  pressure  correlation 

Mass  of  dissolved  solids  equivalent  to  a unit  mass  of  TOC  for  use  in  the 

osmotic  pressure  correlation 

Coefficient  in  mass  transfer  correlation  for  A 

Coefficient  in  mass  transfer  correlation  for  C 


COST  FUNCTION 

The  test  data  basically  consists  of  two  types: 

1.  Permeate  flux  data 

2.  Permeate  concentration  data  (dissolved  solids  and  TOC) 

The  above  values  are  read  directly  or  computed  from  the  experimental  data. 

The  cost  function  employed  is  the  following: 


Table  111.  Comparison  of  Experimental  Data  and  Model  Predictions 


Permeate  Concentrate 


Flow  rate  (irr/hr)  1.02 
Solids  concentration  (gm/m3)  12100.0 
TOC  concentration  (gm/m3)  581.0 


0.36 

104.0 

31.0 


0.66 

18600.0 

881.0 


Permeate 

0.335 

105.0 

30.0 


where : Zi*  WFi*  WDi*  and  Wj...  are  weighting  factors  for,  respectively,  the  point  as  a 
whole,  permeate  flux,  permeate  solids  concentration,  and  permeate  TOC  concentration.  * 

^A’  **C  rePreserT^  permeate  dissolved  solids  and  TOCconcentrations,  respectively.  Subscri; 
e and  m indicate  experimental  data  and  model  predictions,  respectively. 

Table  III  provides  a comparison  between  a typical  data  point  and  model  predictions  for  tl 
same  feed  and  operating  conditions  after  having  determined  values  of  the  model  parameter: 
as  outlined  above.  (Module  pressure  = 54.37  atm,  temperature  = 302.78°K) 


Concentral 

0.61 

18000.0 

850.0 
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A Mathematical  llodel  of  an  Ultrafiltration 
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Abstract 

This  paper  describes  the  development  of  a mathematical 


model  depicting  the  ultrafiltration  unit  employed  in  the  water 
processing  element  under  development  by  the  U.S.  Army  Medical 
Research  and  Development  Comand  for  support  of  field  medical 
units.  Methodology  for  solving  the  model  equations  and  fitting 
the  model  parameters  to  experimental  operating  data  has  been 
developed  and  programs  written  for  a digital  computer.  This 
article  presents  in  detail  the  model  equations,  the  solution 
techniques  employed,  the  parameter  fitting  procedures,  and  some 
typical  results. 


Development  of  Basic  Equations 

The  model  equations  developed  charaterize  the  operation 
of  the  Abcor  HFO  ultrafiltration  unit  shown  in  Fig.l.  Contami- 
nated water  enters  the  tube  at  z* 0 and  flows  in  the  z direc- 
tion, exiting  atz=L,and  is  collected  as  concentrate.  While  in 
the  tube,  some  of  the  water  and  contaminants  pass  through  the 
membrane  attached  to  the  tube's  Inner  surface  and  is  collected 
as  permeate. 


:Ao 


'•Bo 


. 1 

LZ  1 

c 

F(z) 

F(z+az)  } ^ 
» D_ 

V J 

C.(z)  ’ 

C.(z+Az)  I1" 

v 

' C*(Z) 

c'b(^z>|  J 
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LI  ** 
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Fig. 
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1.  Abcor  MFD  Tubular  Ultrafiltration  Unit 
inlet  concentration  of  contaminant^  gm/m3 
inlet  concentration  of  water;  gm/m 


- diameter  of  module;  m 
F„  - Inlet  flow  rate;  m3/hr 


o ■ 2 

J.  - flux  of  contaminant  across  membrane;  <y>/ * -hr 
M a 


flux  of  water  across  membrane;  gm/m  -hr 


L - length  of  module;  m 
z - axial  spacial  coordinate;  m 


for  the  present  work,  the  model  Is  only  concerned  with  the 
suspended  solids  and  the  dissolved  solids  In  the  feed.  Experi- 
mental data  Indicates  that  both  have  an  effect  on  the  flux  of 
water  and  dissolved  solids  through  the  maabrane,  the  permeate 
being  essentially  free  of  any  suspended  solids.  As  presently 
written,  the  model  assumes  that  the  effect  of  suspended  solids 
and  dissolved  solids  on  the  flux  Is  the  same,  which  permits  the 


dissolved  and  suspended  solids  to  be  luaped  Into  a single  eon- 

S.  Ml 


tamlnant,  'A*,  for  the  feed  and  concentrate  streams. ... 

experimental  data  becomes  available,  It  may  be  possible  to  dis- 


tinguish between  the  effects  of  dissolved  solids  and  suspended 
solids  on  the  fluxes. 


A steady  state  material  balance  on  the  total  solids,  compo- 
nent 'A',  for  the  differential  element  in  Fig.  1 yields; 


Input  rate  of  A * Output  rate  of  A 

F(z)  CA( z)  • F(z+  az)CA(z*  bz)  * nDfflJ^(z)AZ 


(1) 


where 


: F(z)  ■ 
C.(z)  - 


fluid  volumetric  flow  ra^e  (m/hr) 


concentration  of  A (gm/m3 ) 


JA(z) 


flux  of  A through  membrane  (gm/hr-m  ) 


As  the  permeate  contains  no  suspended  solids,  J.  must  consist 
entirely  of  dissolved  solids. 


Similarly,  a steady  state  total  material  balance  on  the  dif- 
ferential element  in  Fig.  1 yields: 


Input  rate  * Output  rate 


where: 


F(z)p  ■ F(z+  Az)p  ♦ vDBtJA(z)  ♦ Jg(z) )Az 

.3, 


(2) 


total  mass  density  (gn/m  ) 
flux  ‘ 

(gm/m 


J.(z)  * flux  of  water,  component  B,  through  the  membrane 
' i/mZ-hr) 


In  addition,  a steady  state  material  balance  on  component 
'B'  can  be  written  to  obtain  an  equation  analogous  to  equation 
(1),  However,  this  equation  is  not  necessary  as  it  becomes  a 
dependent  equation  when  equations  (1)  and  (2)  are  used  in  the 
model . 


Equations  (1)  and  (2)  are  readily  converted  to  differential 
equations  by  taking  the  limit  as  Az  approaches  zero.  For  equa- 
tion (1).  the  result  Is: 


ar  if(*)caCz)j  * - -vA<z) 

Equation  (2)  becomes: 

nD_ 


(3) 


-P-  * Vz)1 


(4) 


Each  ultrafiltration  tube  Is  modeled  as  a long  hollow  cy- 
linder with  a semi-permeable  membrane  adhering  to  the  Inner  sur- 
face of  the  tube.  The  membrane  acts  as  an  extra  fine  mesh,  re- 
taining suspended  solids  and  some  very  large  organics,  but  pass- 
ing water  and  such  of  the  dissolved  solids.  During  operation, 
a thin  boundary  layer  Is  asstsaed  to  form  at  the  membrane  through 
triilch  water  must  pass  before  reaching  the  tube  wall.  Detailed 
equations  derived  from  writing  component  and  overall  material 
balances  for  the  bulk  and  boundary  layer  sections  of  the  tube 
form  the  basis  of  the  model . 


The  boundary  conditions  are  as  follows: 


Fo<°>  * Fo 


CA<°)  ' CAo 


where: 


F • feed  rite  to  nodule,  m /hr 
o 


• Inlet  concentration  of  A, 


nU,*U,^nirftf5;  and  0,  a stagnant  boundary  layer  1, 
Imagined  to  exist  near  the  inner  will  of  the  tube  as  shown  In 
Fig.  2.  The  following  fluxes  Into  and  out  of  the  boundary  layer 
are  defined: 


gDA 


»<CA2  * CA3> 

(S) 

y(AR  - Av) 

(*) 

kA(CA2  ‘ CAJ* 

m 
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JDB  “ kB*CBl  ‘ CB2* 


Bulk  Fluid 


cAi,cbrxArxBi 


J0A  I JDB 

-1-1-4 


Boundary  layer 


r r x x 
-A2,V*B2*  A2*  B2 


Permeate 


JA  JB  CA3,CB3’XA3’XB3 


Fig.  2.  Fluxes  into  and  out  of  Boundary  Layer 


CAl  “ concentration  of  A in  bulk  fluid,  yn/m3 
*"A2  " conc*ntration  of  A in  boundary  layer,  gm/m3 


C*3  * concentration  of  A in  permeate,  gm/nr 


Cgj  * concentration  of  B in  bulk  fluid,  gm/m 


Cg2  * concentration  of  B in  boundary  layer,  gm/m 
Cg3  ■ concentration  of  B in  permeate,  gm/m3 


y * permeability  coefficient  for  pure  water,  yn/m- 
hr-atm 

B ■ permeability  coefficient  for  dissolved  solids, 
m/hr 

bP  * pressure  drop  across  membrane,  atm 

bn  - osmotic  pressure  difference  across  membrane  * 


*2  - »j,  atm 


*2  * osmotic  pressure  at  boundary  layer  fluid,  atm 
Vj  * osmotic  pressure  of  permeate,  atm 


k,,k.  • mass  transfer  coefficients  for  solids  and  water, 
respectively,  m/hr 


The  osmotic  pressure  Is  related  to  the  concentration  of 
contaminant,  CA  and  the  absolute  temperature,  T,  by  the  empiri- 
cal relationship 

* * f(CA.T)  * aTCA(1  * KA)2  (9) 


where:  a ■ coefficient,  atm-m /g*.  K 
B * coefficient,  m^/gm 
T • absolute  temperature,  deg  Kelvin 


To  obtain  be  • ».  - w.,  It  is  necessary  to  evaluate  «,  using 
CM  and  *j  uslng'c^j. 


The  flux  J.  due  to  the  bulk  flow  toward  the  boundary  layer 
can  be  obtainedrby  writing  a total  material  balance  around  the 
boundary  layer: 


JF  ' JA  * J»  * JDA  ‘ JDS 


In  addition,  a component  balance  for  component  A con  be  written: 


JF*A1  * JA  + J0A 


where  XM  ■ mess  fraction  of  dissolved  solids  In  the  bulk  str 


A component  balance  can  also  be  written  for  component  I,  but  It 
will  not  be  an  Independent  equation. 


In  addtion,  X.,.  the  mass  fraction  of  dissolved  solids  in 


the  permeate  is  related  to  and  Jg  as  follows: 


This  completed  the  set  of  equations  describing  the  fluxes  at  the 
boundary  layer. 


In  the  ultrafiltration  unit,  the  concentration  of  dissolved 
solids  will  be  relatively  small.  Furthermore,  the  membrane  will 
tend  to  reject  these  components.  For  these  reasons,  the  follow- 
ing assumptions  will  be  made: 


1.  The  membrane  permeability  flux  is  small  compared  to 
JB‘ 


2.  The  diffusional  fluxes  and  JnR  are  small  compared 


3.  The  density  is  that  of  water. 


where:  JM  - ^/fusionaW lux  of  contaminant  out  of  the  boundary 
JoB  " gm/m^hr"41  flUX  °f  "ater  1nt0  the  boundarX  layer. 


Table  I presents  the  flux  equations  with  these  approximations 
Included.  Observe  that  it  is  no  longer  necessary  to  evaluate 

JDB' 


Simplified  Equations  for  Fluxes  at  the  Boundary  Layer 


Transport  Equations 


Osmotic  Pressure  Relationships 


B(CA2  * CA3) 


J.  • y(bP  - be) 


*2  * oTCA2(1  4 *CA2) 


JDA  * kA<CA2  - CA1> 


*3  ‘ oTCA3(1  + 6CA3' 


Material  Balance  at  Boundary  Laver 


JA  4 3 DA 


x • 

XA3  r 


Mass  hractlon/Density/Concentration  Relationship 


XB2  * 1 * XA2 


*83  * 1 * XA3 


C,,  * PX. 


CA3  * **A3 


CB2  * ^62 


CB3  * 0*83 


Solution  Method 

The  solution  method  basically  consists  of  numerically  Inte- 
grating the  differential  equations  (3)  and  (4)  resulting  from 
the  total  and  component  material  balances.  However,  this  Is 
somewhat  complicated  by  the  fact  that  the  fluxes  must  be  deter- 
mined at  each  Increment  by  solving  the  set  of  nonlinear  equa- 
tions summarized  In  Table  I. 


The  computation  of  J,  and  J.  at  every  point  z re- 
quires the  solution  of  tM  set  of® equations  in  Table  I.  which 
entails  two  nested  Iteration  loops. 


For  both  loops,  the  Interval  halving  method  was  found  to 
be  necessary  to  achieve  a reliable  convergence  of  the  loops.  A 
similar  situation  exists  In  the  reverse  osmosis  unit  downstream 
of  the  ultraflltrelon  unit.  The  equations  proved  to  be  highly 
nonlinear,  and  with  varlng  degrees  of  sensitivity.  Consequently 
the  reliability  of  the  Interval  halving  method  had  to  take  pre- 
cedence over  computational  efficiency. 


For  the  first  or  outer  iteration  loop,  a value  of  must 


•\v 
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be  estimated.  It  Is  known  that  CAJ  £ £ o.  which  permits  in- 

terval halving  to  be  utilized. 

However,  knowing  the  value  of  only  permits  to  be 
calculated.  To  proceed  further,  a value  is  assumed  for  C^j, 
which  must  be  in  the  interval  0 < C«,  < C...  This  permits  J.  to 
be  readily  calculated  as  follows: 


(13) 


by  combining  the  following  three  equations: 


If  C^*  is  not  equal  to  C^,  a new  value  for  C^2  must  be  assumed. 

The  Interval  halving  logic  permits  the  interval  of  uncertainty 
for  CA2  to  be  reduced  and  the  calculations  repeated. 

Using  this  procedure,  the  nonlinear  flux  equations  are 
solved  for  each  Integration  step. 

Parameter  Estimation  Considerations 

Values  of  several  key  variables,  such  as  mass  transfer  co- 
efficients and  osmotic  pressure  correlation  coefficients,  must 
be  determined  from  experimental  data.  A Pattern  Search  strat- 
egy (4)  was  used  to  find  a set  of  values  for  these  key  variables 
which  minimized  the  difference  between  model  calculations  and 
experimental  results. 
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Harriott  and  Hamilton  (2)  give  the  following  relationship 
for  tubular  ultrafiltration  modules: 


0.913  H 0.346 

"sh  * °-0096  nr.  Sc 


where:  Ngh  * 


kD 

PP  • Sherwood  number 
°AB 


(14) 


The  Pattern  Search  technique  was  chosen  for  several  rea- 
sons: 

• Simplicity 

• Ease  of  use 

. Relative  insensitivity  to  numerical  error,  i.e., 
no  partial  derivatives  are  evaluated 

• Reasonably  fast  convergence  to  the  solution 
The  method  is  fairly  well  known 

Table  II  lists  the  parameters  whose  values  have  been  obtained 
by  the  parameter  estimation  procedure. 


DV 
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Reynolds  number 


N,  * x— - — * Schmidt  number 
ic  UAB 

where:  k • mass  transfer  coefficient,  m/hr 
v * kinematic  viscosity,  m^/hr 
O^u  - diffusivity,  r.i^/hr 

D * characteristic  length,  m 

This  relationship  was  Incorporated  directly  into  the  model. 
However  in  attempting  to  improve  the  accuracy  of  the  model,  the 
relationship 


UAB 


a* 


A1  > b 


was  used  to  make  the  diffusivity  a function  of  concentration 
where  a,  b,  and  c are  coefficients. 

Now  that  values  are  available  for  and  CgJ.  the  osmotic 
pressures  x^  and  Xj  can  be  computed,  and  then  Jg*  can  be  evalu- 
ated using 


and 


Ax  * x,  - x. 


Jg*  ■ (AP  - Ax) . 


If  the  values  of  Jg  and  Jg*  do  not  agree,  the  value  of  must 

be  changed.  Using  the  interval  halving  technioue,  the  interval 
of  uncertainty  for  Cgg  is  reduced  and  the  computations  repeated. 

Once  a value  of  C.,  is  obtained  that  gives  the  same  value 
for  Jg  and  the  assumption  for  must  be  checked.  Since 

values  are  available  for  JB,  J,,  and  J-, , we  can  compute  C,-*  as 
follows:  B A DA 


„g,  -g,  xM, 

VaI  ’ kA(CA2*  - CA1»  + "<CA2*  ' CA3> 


rewritten  as 


jbxai  * kac, 


‘•A2 


T^TT 


y*Ai  * Ku 


Table  II 
Fit  Parameters 

y Permeability  coefficient  for  pure  water 

B Permeability  coefficient  for  A 

a Coefficient  in  osmotic  pressure  correlation 

B Coefficient  in  osmotic  pressure  correlation 

a,b,c  Coefficients  in  the  diffusivity  relationship 


The  experimental  data  on  which  the  parameters  were  fit 
was  essentially  a “snapshot"  of  a steady-state  condition  for 
the  UF  unit.  The  following  Information  was  available: 

1.  Specifics  of  Module  Configuration  and  Geometry 

a.  Length  (m) 

b.  Inner  diameter  (m) 

c.  Inner  surface  area  of  module  available  for  flux 
(m2) 

2.  Operating  Test  Data 

a.  Inlet  volumetric  flow  rate  (m3/hr) 

b.  Inlet  contaminant  concentration  (yn/m3) 

c.  Inlet  pressure  (atm) 

d.  Outlet  pressure  (atm) 

e.  Temperature  (°K) 

f.  Average  flux  of  contaminant  (gn/n^-hr) 

g.  Average  flux  of  water  (gm/m2-hr) 

Table  III  lists  one  set  of  the  experimental  data  that  was  avail- 
able. 
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Table  111 
Experimental  Data 


Diameter:  0.0254  n 

2 

Inner  Surface  Area:  1.64  m 
Length:  24.4  m 


where:  COST  « cost  function  to  be  minimized 
1 * date  point  Identifier 
n • number  of  data  points 


z,  » weighting  factor  of  point  1 relative  to  other 
1 data  points 


Wci  * weighting  factor  for  flux  for  data  point  1 

F1  (0<wF1<i) 


Wn.  - weighting  factor  for  contaminant  concentrator 

^ ' Cab  <4af  s nninf  < I f\  < ^ 1/ 


u for  data  point  1 (0  < W0<  <,  1) 

F . • experimental  value  of  permeate  flux  for  data 
8,1  point  1 

F^  1 ■ permeate  flux  predicted  by  model  for  data  point  1 
C.  . « experimental  permeate  dissolved  solids  concentra- 
tlon  for  data  point  1 

C,  . • permeate  contaminant  concentration  predicted  by 
*"•'  model  for  data  point  1 


In  addition  to  the  requirement  that  the  weighting  factors  WF( 
and  WD1  must  each  be  less  than  1.0,  they  should  logically  st«i  to 
1.0  for  each  point  1: 


WF1  * M01 


Table  IV  Illustrates  the  quality  of  the  fit  to  the  data 
given  In  Table  III.  The  weighting  factor  z,  was  1.0  for  each 
point,  with  M,.  and  U„4  being  0.8  end  0.2,  Respectively,  for 
each  point.  M 


Table  IV 
Model  Fit 


Permeate  Flux 


ate  Concentration 
/ _ ,—3 1 


Continuing  Work 

Present  efforts  are  being  directed  toward  Incorporating 
TOC  (Total  Organic  Carbon)  relationships  Into  the  model,  but 
work  Is  awaiting  the  availability  of  experimental  data  on  TOC. 

At  that  time,  additional  efforts  will  be  directed  toward  improv- 
ing the  model  fit  beyond  that  Illustrated  In  Table  IV.  The 
ultrafiltration  model  Is  also  being  combined  with  models  for 
reverse  osmosis,  ozonation,  and  hypochlorlnatlon  to  provide  a 
comprehensive  model  for  the  water  processing  unit. 
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Inlet  Concentration 

Permeate  Flux 

Permeate  Concentration 
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